Fast-ion conduction and flexibility and rigidity of solid electrolyte glasses 
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Electrical conductivity of dry, slow cooled (AgP03)i-i:(AgI)i glasses is examined as a function of 
temperature , frequency and glass composition. From these data compositional trends in activation 
energy for conductivity Ea(x), Coulomb energy Ec(x) for Ag^ ion creation, Kohlrausch stretched 
exponent /3(x), low frequency (es(x)) and high-frequency (eoo(x)) permittivity are deduced. All 
parameters except Ec(x) display two compositional thresholds, one near the stress transition, x = 
Xc(l)= 9%, and the other near the rigidity transition, x = Xc(2)= 38% of the alloyed glass network. 
These elastic phase transitions were identified in modulated- DSC, IR reflectance and Raman scatter- 
ing experiments earlier. A self-organized ion hopping model (SIHM) of a parent electrolyte system is 
developed that self-consistently incorporates mechanical constraints due to chemical bonding with 
carrier concentrations and mobility. The model predicts the observed compositional variation of 
(t(x), including the observation of a step-like jump when glasses enter the Intermediate Phase at 
x>Xc(l), and an exponential increase when glasses become flexible at x>Xc(2). Since Ec is found 
to be small compared to network strain energy (Ea), we conclude that free carrier concentrations 
are close to nominal Agl concentrations, and that fast-ion conduction is driven largely by changes 
in carrier mobility induced by an elastic softening of network structure. Variation of the stretched 
exponent /3(x) is square- well like with walls localized near Xc(l) and Xc(2) that essentially coincide 
with those of the Intermediate Phase (IP) (xc(l)<x<Xc(2)), and suggest filamentary (quasi ID) 
conduction in the IP, and conduction with a dimensionality greater than 1 outside the IP. 

PACS numbers: 61.43.Fs-61.20.-x 



I. INTRODUCTION 

Fast-ion conduction in solid electrolyte glasses has 
been investigated^ for over two decades and our ba- 
sic understanding continues to be challenged with new 
findings^-^ in the field. Resolution of these basic 
challenges'^'^ including the role of molecular structure and 
impuritiesii on fast-ion conduction is likely to have a di- 
rect impact on applications of these materials in solid 
state batteries, sensors, and non-volatile memories. 

Electrical conductivity a can be generally be written 

as 

a = Ze^riL (1) 

where the terms Ze, il, til represent respectively the 
charge of the conducting ion, its mobility and concentra- 
tion. Beyond this point of discussion, it is clear that al- 
loying ionic species in a base glass will lead to an increase 
in carrier concentration in a host network, but it will also 
substantially modify the mechanical behaviour by either 
depolymerising or softening the base structure. Solid 
electrolyte additives such as Agl, Ag2S and Ag2Se in 
base chalcogenide glasses are known to possess a low con- 
nectivity and their glass forming tendency derives from 
that feature of structurei^^. Recently Ingram et al. 
have proposedi^ that the activation energy for conduc- 
tion (Ea) could be displayed under the form of a master 



plot with respect to activated volume (V a) as 

Ea = MVa (2) 

where M represents a localized modulus of elasticity 
(strongly influenced by intermolecular forces) with val- 
ues ranging between 1 and 8 GPa depending on the na- 
ture of chemical bonding involved. One may, therefore, 
expect that ionic conduction can be related to an elastic 
softening of a glass network with one functionality (elas- 
ticity) affecting the other (conductivity). Can this view 
be made more quantitative? 

Conductivity thresholds in fast ionic conductors have 
been reported^"* although the relationship with rigidity 
transitions was not made. Micro-segregation effects have 
also been reportedi^ in various modified oxide glasses 
including silicates and supported^® by emerging length 
scales in MD simulated partial pair distribution func- 
tions. These findings are popularized by the channel 
picture for alkali ion-conduction in sodium silicates due 
to Greaves and Ngaii^. In this glass system, the soda 
content where a micro-segregation onsets is nearly the 
same as the stress elastic phase transition (x = 18%) 
determined from calorimetric measurementsl-J. A natu- 
ral question that follows then is, are there links between 
the elastic behaviour of glasses and micro-segregation ef- 
fects? Is this a generic behaviour or is it typical of alkali- 
modified silicates? 

(AgP03)i_a;(AgI)a; glasses have been widely studied 
in the past by several groupsi^^. A perusal of this 
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FIG. 1: (Color online) Plot of electrical conductivity (a) mul- 
tiplied by T as a function of 1000/T at the three indicated 
Agl content x in (AgP03)i_a;-(AgI)a: glasses displaying an 
Arrhenius variation. 

large body of work reveals broad compositional trends 
in thermal, optical and electrical properties of these ma- 
terials, but not without significant variations between 
various groups. In some case the conductivity increases 
monotonically with Agl content displaying either no 
threshold^^ or one threshold near x = 30%21-^. In more 
recent reports^J^i^i^ we have emphasized the need to 
keep samples dry during synthesis and handling, and to 
slow cool them from Tg to assure measuring what we 
believe are the intrinsic physical properties of these ma- 
terials. The thermal, optical, molar volumes and some 
electrical properties were reported^S., and these are found 
to be quite different from those reported earlier in the 
literature^Si^^iSa in that many properties including con- 
ductivity display two striking compositional thresholds. 
The data from these more recent investigations^''^'' ex- 
amined as a function of glass composition are found to 
display striking correlations between various observables. 
For example, these data have permitted connecting as- 
pects of local structure deduced from Raman and IR vi- 
brational spectroscopy to their thermal as well as elec- 
trical transport properties. The elastic nature of back- 
bones derived from Lagrangian constraints, have permit- 
ted treating short-range intermolecular forces in an effec- 
tive way using rigidity theory. 

In the present contribution we present temperature de- 
pendent and frequency dependent conductivity data on 
these glasses as a function of Agl content x. Here we 
address an issue of major importance in the field of solid 
electrolytes, namely, the relationship of the three defined 
elastic phases: flexible, intermediate, stressed rigid to 
the dielectric permittivity derived from the frequency- 
dependent electrical modulus. We find that (a) the di- 
electric permittivity is a maximum in the intermediate 
phase (IP), (b) the Kohlrausch stretched exponent /3 dis- 



plays a minimum there and (c); we also infer there are 
dimensional changes in conduction pathways with chang- 
ing glass composition, and (d) that Coulomb barrier , Ec, 
determined from the high-frequency permittivity makes 
only a small contribution to the Arrhenius activation en- 
ergy for dc conduction, Eyi, compared to network strain 
energy, E^. These findings suggest that free carrier con- 
centrations are nearly equal to one, and that they closely 
follow the Agl concentrations. The observed composi- 
tional trends in dc conductivity, (t(x) appear to derive 
from elastic network softening which controls carrier mo- 
bility. The observed precipitous increase of conductivity 
in the flexible phase may well be a generic feature of solid 
electrolyte glasses. 



II. EXPERIMENTAL 
A. Synthesis of dry AgI-AgP03 glasses 

Bulk (AgP03)i_a;(AgI)a; glasses were synthesized us- 
ing 99.9% Ag3P04 (Alpha Aesar Inc.), 99.5% P2O5 ( Fis- 
cher Scientific Inc.), and 99.99% Agl (Alpha Aesar Inc.), 
and separately AgNOa and NIl4H2(P03)4 with Agl as 
the starting materials. Details of the synthesis are dis- 
cussed elsewhereii. In our experiments, the starting ma- 
terials were handled in a dry nitrogen gas purged glove 
box (Vacuum Atmospheres model IIE-493/MO-5, rela- 
tive humidity ^ 0.20%). The starting materials were 
weighed in the desired proportion and thoroughly mixed 
in an alumina crucible with all handling performed in 
the glove box. Mixtures were then transferred to a box 
furnace held at 125°C in a chemical hood purged by lab- 
oratory air, and heated at IOO°C/hr to 700°C. Melts 
were equilibrated overnight and then quenched over steel 
plates. Calorimetric experiments also revealed^^ that 
stress frozen upon melt quenching can be relieved by 
cycling samples through the glass transition tempera- 
ture (Tg). Conductivity measurements were undertaken 
only on samples cycled through Tg and slow cooled at 
3°C/min. to room temperature. 



B. AC conductivity measurements as a function of 
temperature and frequency 

Glass sample disks about lOmm in diameter and 2 
mm thick were synthesized by pouring melts into spe- 
cial troughs. Pellets were then thermally relaxed by 
cycling through Tg. Platelets then were polished, and 
the Pt electrodes deposited. The complex impedance 
Z*(a;)=Z'((jj)-|-iZ" (w) of the specimen was measured by 
a Solartron SI 1260 impedance analyzer over a fre- 
quency range of 1 Hz-1 MHz (19 points per measure- 
ments). The temperature was raised from 150 K to Tg 
at a rate of l°C/min. All data were acquired auto- 
matically at 2 min intervals. The complex conductiv- 
ity a* {uj)=a\uj)+ia" (lo) is deduced from the complex 
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FIG. 2: (Color online) Variations in activation energy Ea(x) 
(blue) and preexponential factor (To(x) (red) as a function of 
Agl content x in percent in (AgP03)i-i:-(AgI)i, glasses. 

impedance Z*, the sample thickness t and the area S 
of the surface covered by platinum using the formula : 
a*(c.)=(t/S)(l/Z*). 



III. EXPERIMENTAL RESULTS 

A. Temperature dependence of conductivity 

The temperature dependence of electrical conductivity 
for ionic conduction in glasses usually displays an Arrhe- 
nius dependence 

CT = (To exp[-£'^/fc_BT] (3) 

We have obtained the activation energy FiAix), by per- 
forming temperature dependent measurements in the 
40°C<T<Tg range. An illustrative example of some of 
these results, at three glass compositions, is reproduced 
in Fig 1. From the temperature dependence, we have 
obtained the preexponential factor (Tq a-nd the activation 
energy as a function of Agl content, and these data 
are reproduced in Fig. 2. One finds that E^(x) steadily 
decreases as x increases, particularly for x > 40% while 
CTo drops rapidly with x at low x (<10%) and then levels 
off in the intermediate x region. 

Interestingly, we find cr(x) displays three regimes of 
conduction whatever the temperature (Fig. 3a). At low 
x (< 9%) cr(x) is found to be nearly independent of x. 
The conductivity then builds up at higher x (> 10%) 
with a distinct step near 9%, and at x > 37.8% a power- 
law variation of (t(x) sets in. The location of the two 
electrical conductivity thresholds at x = 9% and x — 
38% coincide with steps in the non-reversing enthalpy at 
Tg that define the boundary of the Intermediate Phase 
(IP) (Figure 3b). These data will be discussed in section 
IV. 
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FIG. 3: (Color online) a: Variations in DC conductivity cr(x) 
at three temperatures, -50°C, 50°C and 100°C, and b: Non- 
reversing heat flow AH„r at Tg in (AgP03)i_2,-(AgI)a; glasses 
as a function of Agl content in mole %. 



B. Frequency dependence of AC conductivity 

The frequency dependence of AC conductivity, <t{uj) at 
three representative glass compositions is summarized in 
Figure 4. These data show a saturation of cr(w) at low 
frequencies ([x)/27r < 100 Hz), and an increase of about 
2 orders of magnitude in the 10'^ < cu < 10^ range. The 
low frequency saturation value of a is found to increase 
steadily as the Agl content of glasses x increases from 
10% to 40%, a behaviour that is consistent with data 
in Fig. 3. The present results on the frequency depen- 
dence (j{uj) are similar to earlier reports of Sidebottom^ 
and separately by Stanguenec and EUiott^^. For a more 
quantitative analysis of these a{uj) data, we have deduced 
the frequency dependence of the dielectric permittivity, 
eoo(w) from the measured AC conductivity. 



C. Electrical modulus analysis 

The behaviour of the complex conductivity 
a* {uj)=a\u!)+ia" (uj) can be analyzed using the complex 
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FIG. 4: Frequency dependence of electrical conductivity in 
(AgP03)i_i-(AgI)a; glasses at three indicated Agl concentra- 
tions X. Measurements were taken at 250K. 



permittivity £*(w)=e'(aj)+i£"(a;) or the complex electri- 
cal modulus M*(tj)=M'(w)-HiM"(tj) which are related 
to the complex conductivity via: 



and 



M*iuj) 



ilusq 



(4) 



(5) 



Here Sq designates the permittivity of free space of 
8.85 X 10~^^ Farad/meter. From the measured value of 
cr" , using equation (4) we could deduce the real part of 
the permittivity e'. These data are presented on a log- 
log plot for a glass at a composition x — 25% in Figure 
5. Permittivity results generally from mobile ions, but 
the sharp increase in e' at low frequencies, particularly 
at f = a)/27r < 10 Hz, is due to electrode polarization ef- 
fects as Ag+ ions are deposited at the anode at very low 
frequencies'^^ . The behaviour is observed for all samples 
but the magnitude of £s could be unambiguously deter- 
mined. The saturation value of the permittivity at low 
frequency (Es) and at high frequency (Eoo) were thus ob- 
tained at all other glass compositions, and these data ap- 
pear in the insert of Fig 5. Noteworthy in these trends of 
Eg (x) and £00 (x) is the fact that both show a broad maxi- 
mum that onsets near x = xi ~ 10% and ends near x = X2 
~ 35%, the two threshold compositions observed in the 
non-reversing heat flow at Tg (Fig. 3b) and electrical con- 
ductivity (Fig. 3a) that essentially define the IP region^ of 
these glasses. At the glass composition x = 20% , located 
near the middle of the IP, these parameters acquire rather 
large magnitudes: £^=620 and £oo=250. The magnitude 
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FIG. 5: (Color online) Log-log plots of the real part e' of per- 
mittivity vs frequency for (AgP03)i-a:-(AgI)a; glasses at x = 
25% measured at 250 K. The broken horizontal lines corre- 
spond to the high and low-frequency limits of the permittivity 
£00 and Es- The inset shows variation of £00 (x) and es(x) with 
Agl content. Vertical broken lines define the location of the 
stress and rigidity transitions at Xc(l) = 9% and Xc(2) = 38% 
as determined from calorimetry. 



of these permittivity found in the IP of the present elec- 
trolyte are truly exceptional- they are high compared to 
those found in other systems such as the alkali silicates 
(£oo=llr^), and highlight another remarkable property 
of sohd electrolytes bearing Ag salts. 

We now follow the approach developed by Ngai and 
co-workers2i and use the electrical modulus formalism to 
obtain stretched exponents /3(x) from our a{uj) data. Our 
starting point is to write the electric field E(t) under the 
constraint of a constant electric displacement as: 



E{t) = E{Q)^t) 



(6) 



and define E(0) as the initial electric field and $(t) the re- 
laxation function of the system. In the frequency domain, 
the electrical relaxation is appropriately represented by 
the electric modulus: 
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(7) 



and the high-frequency (£oo)and low-frequency limits [Es) 
of the real part of the permittivity (£ ) then can be writ- 
ten as : 



< > / < r >^ 
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FIG. 6: Variations in the Kohlrausch stretched exponent /3(x) 
as a function of the Agl concentration x in (AgPOa)!-!- 
(Agl)^; glasses. 



where: 



/•CXD 

Jo 



(9) 



We use the standard relaxation function based on a 
stretched exponentiah 



(10) 



with representing the Kohlrausch- WiUiams- Watts expo- 
nent. The exponent provides a good measure of the cor- 
relations of atomic motions in a relaxing network'^^. And 
the absence of cooperative behaviour will lead inevitably 
to a typical Debye relaxation with (3 —1, a behaviour that 
is expected in dilute limit where concentrations of diffus- 
ing cations is low. One finally obtains from equations 
(8)-(10): 



''[r(i//?)]2^= 



(11) 



where F represents the Gamma function. 

The data of Fig. 5 gives the glass composition vari- 
ation of the permittivity, Ssijx) and eoo(x). Using these 
permittivity parameters and the relation between them 
(equation (11)), we have extracted compositional vari- 
ation of the Kohlrausch stretched exponent /3(x) These 
data appear in Figure 6. We find (3 ~ 0.65 at low x 
(0%) and at high x (50%), but at compositions in be- 
tween, particularly in the IP, (3 is found to decrease (/? ~ 
0.53) and to display a broad square-well like minimum, a 
trend reminiscent of the reversibility window (Fig. 3(b)). 
At low carrier concentrations, x ~ 5% for example, we 
are in the dilute limit and carriers can be expected to be 
randomly distributed displaying little or no correlations 
in their motion. Yet our data shows that p does not 
acquire a value of 1, as found in the sodium silicatea^^. 



This is a curious result and it possibly reflects the fact 
that in the present system the diffusing cation of interest 
Ag+ also happen to be part of the base glass network. 



IV. DISCUSSION 

A. Three regimes of fast-ion conduction and elastic 
phases of AgPOs-AgI glasses 

Perhaps the central finding of this work is that fast- 
ion conduction in the present solid electrolyte glass dis- 
plays three distinct regimes of behaviour, and that these 
regimes coincide with the three known elastic phases of 
the backbone. The result emerges from the compositional 
trends in the electrical and thermal data presented ear- 
lier in Figure 3. The observation of a reversibility win- 
dow fixes the three elastic phases of the present glasses^ 
with glass compositions at x < xi — 9% belonging to the 
stressed-rigid phase, those at x > X2 = 38% to the flexible 
phase, while those in between (xi < x < X2) belonging 
to the IP. The two thresholds, xi and X2, we identify re- 
spectively with the stress and the rigidity transitions. As 
glasses become flexible at x > 38%, (t(x) increases loga- 
rithmically, a result that underscores the crucial role of 
network flexibility in opening a doorway for Ag+ ions to 
migrate. In the flexible phase, network backbones can be 
more easily elastically deformed (less cost in strain en- 
ergy), and ions can more easily diffuse and contribute to 
conduction. 

A closer inspection of Fig. 3a also shows that cr(x) is 
nearly independent of x in the stressed-rigid glasses ( x < 
xi), but it varies discontinuously displaying a jump near 
the stress transition, xi ~ 9%. The jump in cr(x) de- 
creases as temperature is increased from -50°C to 100°C. 
At low temperature (-50°C), we note that the observed 
jump is rather large; a increases by a factor 2.5 between 
the composition x =10% and x = 10.5%. The discon- 
tinuity in cr(x) at the stress transition is reminiscent of 
a similar behaviour encountered in the variation of the 
Raman mode frequency (v) of corner-sharing tetrahedra 
in binary Ge^^Si-^: glasses at the stress transition^. The 
frequency jump (Az/) across the threshold also increases 
as T is lowered. Since the stress transition is found to be 
a flrst order transition, it suggests that the first derivative 
of the free energy must be temperature dependent. 



B. Self-organized ion hopping model (SIHM) of 
solid electrolytes 

To gain insights into the origin of the observed electri- 
cal conductivity thresholds (Fig. 3), we use Size Increas- 
ing Cluster Approximation^^^ (SICA) to identify the 
two elastic thresholds in a parent solid electrolyte glass 
system. For this purpose we use the alkali oxide modi- 
fied Group IV oxide, (l-x)Si02-xM20, with M=Li,Na,K, 
and start with N tetrahedra (e.g. Si04/2 and MSi05/2 



(sharing one non-bridging or terminal oxygen) with re- 
spective probabiUties (1-p) and p=2x/(l-x) at the first 
step, 1 = 1, in the agglomeration process. Our goal is 
to compute the probability of finding specific clusters 
in the three phases of interest (fiexible, intermediate, 
stressed rigid). We note that a Si04/2 structural unit is 
stressed rigid (nc=3.67 per atom) while a MSi05/2 unit 
is flexible (nc=2.56 per atom). Here Uc represents the 
count of bond-stretching and bond-bending constraints 
per atom. Henceforth, we denote these local structures 
as St (stressed-rigid) and Fl (flexible) units respectively. 
Starting from these local units as building blocks, we cal- 
culate all possible structural arrangements to obtain ag- 
glomerated clusters containing two (step 1=2), and then 
three building blocks (1=3), etc. and their corresponding 
population probabilities pi. In calculating the probabil- 
ities of the agglomerated clusters we fold in their me- 
chanical energies. Details of the method and application 
appear elsewhere^S,. At each agglomeration step 1, we 
compute the floppy mode count f*^'^. At step 1=2, the 
count of floppy modes f'^' , 

= 3 - = 3 - ^^'^-^^.St^cmPjk ^^2^ 

l^j,k=Fl,St ^jkPjk 

where Tic{jk) £^nd Nj^, are respectively the number of con- 
straints and the number of atoms found in a cluster with 
probability pjk- An IP is obtained if self-organization is 
achieved??"-"^ , i.e. if upon decreasing the content of mod- 
ifler atoms (x), stress- free cyclic (ring) structures form. 
In these calculations we deflne a parameter i] which mea- 
sures the fraction of atoms in ring structures to those 
in dendritric ones. Thus, 77=0 corresponds to the case 
of networks composed of only dendritic structures (no 
rings), while ??=!, corresponds to networks with only ring 
structures. Our calculations reveal that IP-widths are di- 
rectly related to the parameter 77 i.e, larger 77 the greater 
the IP-width. In practice, one starts with a flexible net- 
work containing a high concentration (x) of modifier ions. 
One then decreases the modifier concentration, and can 
expect more St species to form in a cluster that still 
possesses a finite concentration of floppy modes f'^' >0. 
At a certain composition x=Xr, there will be enough St 
species formed to drive clusters to become rigid when 
f*^^^ vanishes, thus deflning the rigidity transition. And if 
we continue to reduce the modifier concentration beyond 
the rigidity transition, one can expect clusters that are 
almost stress- free by balancing St and Fl units. Such a 
selection rule of self-organization will hold up to a certain 
point in composition x=Xs , the stress transition, i.e., at 
X < Xg dendritic or St-rich structures will predominate 
and stress will percolate across clusters. In this SICA ap- 
proach, the IP width is then defined by the composition 
interval, x^-Xg , which may be compared to the experi- 
mentally established interval between the observed stress 
and rigidity transitions, i.e., Xc(2) - Xc(l) = 28% (Fig 3b). 
In what follows, we fix the value 77 in order to work with 
a flxed IP- width. 

Our SIHM theory of ion-conduction builds on an ap- 
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FIG. 7: Structural sites found in the model of (l-x)Si02- 
XM2O: a) a stressed rigid St-St connection (nc=3.67) rep- 
resentative of the base glass Si02, b) an isostatically rigid 
MSi05/2-Si02 connection (St-Fl with nc=3.0) containing one 
M cation, c) a flexible MSi05/2-MSi05/2 connection (Fl-Fl 
with nc=2.56) having two M cations. Both b) and c) can 
give rise to vacant sites such as found on the MSi05/2-Si0^y2 
cluster shown in d) which derives from cluster c). 

proach initially devised for the conductivity of alkali bo- 
rate glasses^^i^^ and identifies hopping sites for conduct- 
ing cations from a structural model, which in our case 
was used to describe application of SICA to (l-x)Si02- 
XM2O glasses. Both (present work and ref^) models in- 
corporate the idea first invoked by Anderson and Stuart^i 
that activation energy for ionic conduction (E^) is made 
of two terms- an electrostatic energy (Ec) to create an 
ion and which determines the free carrier concentration 
n^ , and a network stress energy (Eg) that facilitates mi- 
gration of ions and determines carrier mobility /z. We, 
therefore, concentrate our efforts on estimating the en- 
ergies, Ec and Es, as they can be directly related to the 
statistics of clusters and the enumeration of constraints 
via statistical mechanics averages. 

We begin by modelling the parent solid electrolyte, (1- 
x)Si02-xM20 , where M represents an alkali ion, and 
evaluate the free carrier concentration n^. We will have 
no ions ( n = 0) available if we consider a pair of 
Si tetrahedra, each having one non-bridging oxygen near 
neighbour attached to M+ ions, i.e., a pair of MSi05/2 
tetrahedra (Fl-Fl), since there are no vacancies. How- 
ever, we will have n = 1 and 2 vacancies in the case 
of local structure pairings of the type, MSi05/2-Si0^^2 
, and Si05/2-SiOg^2 i contributing respectively 1 and 2 
M"'' carriers for ion transport (Figure 7). And by know- 
ing the probability for each of these parings pjk, one can 
calculate a mean Coulombic energy over all possible pair- 
ings between Fl and St local units to obtain the carrier 
concentration, 

n^^^ ^2xeyii,[- < Ec> /ksT] (13) 

with the mean Coulomb energy, < Ec > deflned as fol- 
lows, 

<Ec>=^Y.Yl '^^-Pi'^ eM-nEc/keT] (14) 

j,k n 



where Ec is the Coulombic energy to extract a cation 
M+ from an anionic site, and acts as a free parameter in 
the theory. 2, normahses the free carrier concentration 
by requiring n^^ = at T=0, and n^ — 2x at infinite 
temperature, i.e. the nominal carrier concentration. As 
the probability pj^ of pairings depends on the nature of 

the elastic phase, the free carrier concentration n}^ can 
be expected to change as we go across the three elastic 
phases (flexible, intermediate, stressed). 

Once a carrier is available, it will hop between vacant 
sites, and the general form for hopping rates is given by 

= a;^ exp[-i?s/A;BT] (15) 

where is the attempt frequency, and we assume it 
to be constant as only one type of local species (Fl, e.g. 
MSi05/2) is involved in ionic conduction here. In fact, 
given the compositional region where the theory applies 
and where rigidity transitions are usually found^^, one 
does not expect to observe other Fl local species contain- 
ing more than one cation (e.g. M2Si03 corresponding to 
the metasilicate composition). In general, the attempt 
frequency will depend on the local environment of the 
hopping cations'^®. The strain or migration energy Es 
for a hop is roughly the energy required to locally de- 
form the network between a cation site and an available 
vacant site. It should therefore depend on the floppy 
mode energy. We assume it to not depend on the process 
(i,j), which is equivalent to saying that we neglect the 
Coulomb repulsion from vacant sites. Note that in the 
stressed- rigid phase, there are only few hopping events 
possible since the network is made up largely of stressed- 
rigid St-St and isostatically rigid St-Fl pairs. 

In the flexible phase (at x > x^ ) floppy modes pro- 
liferate as the count f^^^ (equation (12)) becomes non- 
zero. In this phase, one can expect local deformations of 
the network to be facile, thus increasing hopping rates, 
and reduce the energy required to create doorways be- 
tween two vacant sites. We write the strain energy as: 
g/iex^gstress_^f(2) ^j^gj.g ^f(2) ^^^^ ^ typical floppy 

mode energy given by experiments^^. In the flexible 
phase, Es is reduced in relation to the stressed-rigid phase 
by the quantity Af'^^^ One is then able to write a con- 
ductivity of the form: 

G = T^l\>P-'> = nf cxp[-S,„/A:bT] (16) 

Results of the self-organized ion-hopping model (SIHM) 
for the parent electrolyte system described above are 
summarized in Figure 8. In this figure we plot the varia- 
tion of conductivity, cr(x) given by equ. (16), as a func- 
tion of glass composition (x) at several temperatures T 
but at a fixed Coulomb energy Ec. This model calcu- 
lation shows that the elastic nature (stressed rigid, in- 
termediate, flexible) of networks alters ion-transport in 
profound ways. Three prominent features become appar- 
ent; (i) a minuscule conductivity change is predicted in 
the stressed-rigid phase, (ii) a step-like increase in (t(x) 
is anticipated at the stress phase boundary x = x^ , and 
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FIG. 8: (Color online) Predictions of the SIHM model for 
variations in ionic conductivity cr(x) for the parent solid elec- 
trolyte, (l-x)Si02-xM20, where M = Alkali ion, using Ec = 
0.1 eV, and at 4 temperatures, from bottom to top: T=300 
K, 500 K, 700 K and 900 K. We have used a medium range 
order (ring) fraction 77=0.6 which defines the IP (see text for 
details) . Curves have been shifted along the ordinate for pre- 
sentation. The broken line corresponds to the conductivity 
when the floppy mode contribution is neglected (see text for 
details) . 



finally (iii) a logarithmic variation of (t(x) is predicted 
in the fiexible phase ( x > x^). These features of our 
SIHM model are in rather striking accord with the ob- 
served variation of conductivity in the present AgPOs- 
Agl glasses (Figure 3). 

The SIHM model predictions provide crucial insights 
into origin of ion-transport in the three elastic regimes. 
In the stressed-rigid phase, strain energy is high and leads 
to a low number of hopping possibilities, since hops are 
allowed only between select number of Fl-St pairs. These 
restrictions on hopping combined with a low free carrier 
concentration results in a weak dependence of (t(x) on 
modifier concentration x. The order of magnitude of 
a normalized to (t(x=0), is largely determined by the 
magnitude of Ec, i.e. the interaction energy between the 
cation M and its corresponding anionic site. 

In the IP, new cation pairs (Fl-Fl, and thus hopping 
possibilities) are populated, and lead to a mild increase of 
a . The step-like increase in conductivity at the stress- 
phase boundary , x = Xg (Fig 7) is a feature that is 
nicely observed in our experiments (Fig 3a) , and provides 
an internal consistency to the description advanced here. 
Theory tells us that the jump in a at stress transition 
depends on T or inverse Ec at fixed T, (see equ. (14)). 
The location of the stress transition is a network prop- 
erty and it depends on the IP, i.e., aspects of medium 
range structure. As in other quantities observed^, the 
first order jump in (t(x) near x = Xs is a manifestation 
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of deep structural changes in a glass network which con- 
tributes to a cusp in the configurational entropy'^'^. And 
as a final comment on the subject, we recognize that at 
a fixed temperature, large values for the Coulomb energy 
Ec will lead, in general, to large jumps in conductivity 
at the stress transition. Here theory predicts that jumps 
in conductivity near the first order stress transition (x = 
Xs) will be more conspicuous with heavier cations (such 
as Ag or K) than with the lighter ones (Li). 

Glasses become flexible at x > x^ as floppy modes 
start to proliferate and decrease the strain energy bar- 
rier. Electrical conductivity displays a second thresh- 
old at the rigidity transition the slope da/dx 
changes abruptly at the transition, and the variation of 
ct(x) becomes exponential above the threshold x^ . Net- 
work flexibility clearly promotes ionic conductivity. In 
the flexible phase an increase in the available degrees 
of freedom appears to facilitate local network deforma- 
tions, thus creating pathways for conduction. Further- 
more, the present results also suggest that the increase 
in the number of possible hopping processes in the IP 
(region between Xj and x^ in Fig. 7) is small and tends 
to produce saturation upon increasing modifler content 
(broken line in Fig. 7). This is to be compared to the 
dramatic increase of a that mostly arises in the flexible 
phase as floppy modes proliferate, f'-^^ ^ , and E^ de- 
creases sharply (Fig 9). 

C. Coulomb energy (E^) and network stress energy 
(E,) in AgP03-AgI 

1. Coulomb energy Ec and Ag^ ion mobilities 

The Coulomb energy Ec represents the energy to cre- 
ate a Ag^ ion that would contribute to ionic conduction 
in the AgPOa-AgI electrolytes glass of interest. In its 
simplest form^^"*,"*" it can be written as 



47r£ooeo RAg-X 

where KAg-x and £00 represent respectively an average 
Ag-X ( X = I, O) nearest-neighbour bond length and 
the high frequency permittivity. From the derived fre- 
quency dependence of the permittivity using the electric 
Modulus formalism, earlier in section IV A, we deduced 
the variation of ecxD(x) in the inset of Fig. 5. The quan- 
tity relates to the dielectric response in the immediate 
vicinity of a silver cation, and permits one to estimate 
the Coulombic barrier needed to create a mobile Ag+ 
cation. In the estimate of Ec we use an average of Ag-I, 
and Ag-0 bond length of 2.43 A and 2.85 A for Raq-x 
in equation (16), and note that the KAg-i bond length 
remains largely constant with glass composition^i. Fur- 
thermore, corrections due to many- body- effects aris- 
ing from Ag^-Ag^ interactions can also be taken into 
account in equation (16) as suggested elsewhere^-. The 
variations of Ec(x) thus largely results from that of £00 (x) 




)l I I ^iJ I I I I I I I I I i_D I I I I L 

10 20 30 40 50 

X (Agl in %) 



FIG. 9: Variations in the activation energy (same as figure 
2), Coulomb energy Ec , strain energy Es — Ea-Ec and equiv- 
alent conductance at 245 K determined from the conductivity 
data of Fig. 3 for (AgP03)i_2!-(AgI)a: glasses. 



(Figure 5 inset), and these are reflected in the plot of Fig- 
ure 9. We flnd that Ec varies between 0.05 eV to 0.10 
eV across the broad range of compositions in the present 
glasses. 

How do the present estimates of Ec compare to earlier 
ones in other systems ? Other groups report similar find- 
ings; in (Na20):E(Si02)i-x glasses. Greaves and Ngaii^ 
estimate Ec near O.leV to 0.2 eV, and the activation en- 
ergies Ea (— 0.6 eV). In (K20)2:(Ge02)i-a; glasses, Jain 
et a\M estimate Ec to be ~ 0.5 eV, typically 30% of E^ 
(~ 1.5 eV) at x < 0.10. Reverse Monte Carlo Models 
of (AgP03)2,(AgI)i_a; glasses by Wicks et al. have also 
emphasized Ec to be small^, and Eg to be the major 
contributor to E^. A similar conclusion is reached by 
Adams and Swenson'*'^, who explicitly reveal that strain 
energy dominates ionic conductivity. The present find- 
ings on Ec highlight once again that the low value of Ec 
appears to be generic feature of solid electrolytes. 

109 NMR studies^ of the base AgPOs glass show 
the NMR resonance to be rather broad, fully consistent 
with its low conductivity (cr ~ 10^^ .cm.~^ , see Fig 
3), in which the Ag"*" cations serve as compensating cen- 
ters. Furthermore, the low value for Ec nearly 10% of the 
total activation energy E^ found in the present AgPOs- 
Agl glassy alloys (Fig. 9) suggests that the number of free 
Ag+ carriers varies linearly with the macroscopic Agl 
content 'x' of the alloys. For this reason, the equivalent 
conductance cr/Zex, i.e. the change in conductivity per 
mole % of Agl, serves a good measure of ionic mobility 
(equation 1). Results on variation of the equivalent con- 
ductance are plotted in Fig. 10, right axis. Here one can 
see that while the nature of the elastic phases weakly af- 
fects Ec, and thus the free carrier density, they induce 
profound variations in the equivalent conductance, and 



thus carrier mobility. In particular, the two thresholds 
in equivalent conductance coincide with those noted ear- 
lier in the variations of Aif„r(x) and a{x) (Fig. 3), which 
serve as benchmarks of the three elastic phases of the 
present glasses. In the flexible phase ( x > 38%) , we 
observe a striking increase by a factor of nearly 20! in 
the 40% < X < 53% range. 

i09^g NMJi studies^" have shown that the number of 
mobile carriers is nearly equal to the number of silver 
ions introduced by Agl. Using the simple relationship be- 
tween a and mobility (Equation (1)) Mustarelli et al.— 
estimated the variation in mobility /i(x) of carriers. Their 
data are also reproduced in Fig. 10 as the red filled cir- 
cles. The corresponding equivalent conductance data on 
our samples are plotted in Fig. 10 as the open square 
data points. It is clear that the mobility in the present 
samples are at least an order of magnitude lower than in 
the samples of Mustarelli et alj^. 

We believe these differences in /z(x) between the two 
sets of samples derive from sample makeup. The glass 
transition temperatures Tg(x) are a sensitive measure of 
bonded water in the base glass, AgPOs, as described 
earlier—. Presence of bonded water in the base glass 
depolymerizes the P-O-P chain network and softens the 
network. The samples of Mustarelli et al,^^ show their 
base glass (AgPOg) T^ = 189°C , about 60°C lower than 
the Tg of our driest base glass of 254°C. Alloying Agl in 
the base glass will soften it further. Furthermore, in our 
samples we do not observe evidence of /3-AgI segregation 
in the 50% < x < 80% composition range as noted by 
Mustarelli et ali^ in their samples. The order of magni- 
tude larger mobilities observed by Mustarelli et al., most 
likely, derive from the flexible nature of glass composi- 
tions promoted by the presence of additional bonded wa- 
ter. 

In our samples, we found^ the IP composition range 
to be rather sensitive to water content of samples. For 
example, the IP-width nearly halved (28% to 15%) and 
the IP centroid moved up in Agl concentration from 24% 
to 28% in going from samples of set A (driest) to those 
of set B (drier). In samples of set B, the T^ of the base 
glass was 220°C, about 34°C lower than in samples of set 
A (Tg = 254°C). And possibly, the solitary conductiv- 
ity (or mobility) threshold observed near x = 30% by 
Mustarelli et al.— (Fig. 10) and separately by Man- 
gion and Johari^^ in their samples may well represent 
the extremal case of a completely collapsed IP at its cen- 
troid near 30% in their samples giving rise to only one 
elastic threshold with the stress- and rigidity- transitions 
merging. These data highlight the importance of sample 
synthesis to establish the intrinsic electrical and elastic 
behavior of these glasses. The role of water traces in 
collapsing the IP in the present solid electrolyte glass is 
reminiscent of a collapse^ of the IP observed in Raman 
scattering of a traditional covalent Ge^jSei-a; glass, mea- 
sured as a function of laser power-density used to excite 
the scattering. In both instances, the self-organized state 
with a characteristic intermediate range order formed, is 
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FIG. 10: (Color online) Equivalent conductance a/Zex (open 
squares, broken line, same as Fig. 7) and silver mobility fj, 
(red circles, right axis) determined from a free carrier rate 
NMR estimation (Mustarelli et al. Ref.^). 



irrevocably changed to a random network structure, in 
one case by OH dangling ends replacing bridging O sites 
and splicing the network, and in the other case by rapid 
switching of lone pair bearing Se centered covalent bonds 
by the action of near band gap light serving as an optical 
pump. 



2. Network Stress energy Es 

Estimates in Ec(x) permit obtaining variations in the 
network stress energy E^ (=Eyi-Ec) from the measured 
activation energies E^(x) (Figs. 1 and 2). Variations in 
the activation energies Eyi(x), and Es(x) are plotted in 
Figure 8 for the reader's convenience. A perusal of these 
data clearly reveals that the drastic increase of electrical 
conductivity in the 40% < x < 60% range is tied to a 
two-fold reduction in Es(x) (Fig. 9). Thus, even though 
the concentration of the electrolyte salt additive (Agl) 
increases linearly with x, changes in conductivity are not 
only non-linear in x, but display discontinuities at the 
stress (xc(l) = 9%) and the rigidity (xc(2)= 38%) elastic 
phase boundaries as discussed earlier. 

The central finding to emerge from an analysis of our 
electrical conductivity results is that fast-ion conduc- 
tion in dry AgPOa-AgI glasses is largely controlled by 
the elastic nature of their backbones. Recent work on 
alkali-earth silicates also shows evidence of conductivi- 
ties increasing precipitously once glasses become elasti- 
cally flexible at a threshold additive concentration^!^. 
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D. Kohlrausch stretched exponent and network 
dimensionality in AgPOs-AgI glasses 

The Kohlrausch exponent usually reflects the degree of 
cooperativity of mobile ions in the glass, the smaller the 
P the larger the collective behaviour of cation motion. 
From what is seen in Fig. 6, one can state that collective 
behaviour is enhanced in the intermediate phase. 

We note that variation of the Kohlrausch exponent 
/3(x) displays a trend that is similar to that of the non- 
reversing heat flow (Fig. 3b) , with values of /3 of about 
0.66 in the stressed-rigid and flexible phases, and a value 
of about 0.52 in the IP. These /3(x) data contrast with 
previous findings of the exponent in the Ge2;Asa;Sei_2a; 
ternary glasses measured in flexure measurements where 
one found''® that /3 steadily increases in going from the 
flexible to the stressed-rigid phases. In the latter it con- 
verged to a value of 0.60 for a network mean coordina- 
tion number larger than 2.4 corresponding to intermedi- 
ate and stressed-rigid phase glasses. 

The stretched exponent /? approaches 1 at high temper- 
atures (T > Tg) , and reduces as T decreases to Tg, i.e., 
/3(T^ Tg) < 1. The reduction of (3 usually reflects con- 
traction of conflguration space in the supercooled liquid 
and its eventual stabilization^ due to structural arrest 
near T=Tg. On this basis at a flxed network dimension- 
ality, one should expect that an increase in the number of 
Lagrangian constraints (or decrease of Agl content) will 
contract the configuration space, leading to lower values 
for p. Fig. 6 obviously does not follow this anticipated 
behaviour as the dimensionality is changing with mod- 
ifier content. Models of traps^^^ for tracer diffusion 
in a d-dimensional lattice show that (3 is related to the 
dimensionality d via: 



f3 = 



d + 2 



(18) 



which is satisfied for various molecular supercooled liq- 
uids and glasses^ having either a 3D (/3 — 3/5) or a 
2D relaxation (/3 — 1/2). However, in the case where 
an internal structural dynamics takes place, the dimen- 
sionality d of networks must be replaced by an effective 
one, deff, involved in relaxation. This happens when 
Coulomb forces are present as the motion of the carriers 
is hindered by the tendency towards local charge neutral- 
ity, determining an effective dimensionality dg// of the 
configuration space in which relaxation takes place^^. In 
this case, the effective dimensionality becomes. 



d = d, 



eff 



(19) 



where is the number of degrees of freedom in labo- 
ratory space and Nc is the number of mechanical con- 
straints, estimated from short-range interactions. In the 
IP where iVc ~ 3, /3 = 1/2 (Fig. 6), equation (17) implies 
d = 2, resulting into a d^ff = 1 from equation (18). Out- 
side the IP, P = 3/5 (Fig. 6), equation (17) implies d — 



3, resulting into a dg// = 1.5 from equation (18). Thus, 
presence of Coulomb interactions in AgPOa-AgI glasses, 
suggests filamentary (dg// = 1) diffusion of Ag+ carri- 
ers in the IP, giving rise to percolative pathways in the 
self-organized structure. This conclusion is supported by 
numerical simulations^, which show that isolated sub- 
diffusive regions exist at low Agl concentration and these 
percolate as filaments at increased Agl concentrations in 
the 20% < x < 30% range. 



V. CONCLUDING REMARKS 

Dry and slow cooled (l-x)AgP03-xAgI glasses synthe- 
sized over a broad range of compositions, < x < 53% 
, are investigated in ac conductivity measurements as 
a function of composition, temperature and frequency. 
From these data compositional trends in activation en- 
ergy for conductivity E^(x), Coulomb energy Ec(x) for 
Ag"*" ion creation, Kohlrausch stretched exponent /3(x), 
low frequency permittivity Es (x) and high-frequency per- 
mitivity eoo(x) are deduced. By combining rigidity the- 
ory with SICA and statistical mechanics, we have de- 
rived the analytic properties of the conductivity response 
function in a Ag based solid electrolyte. The topological 
model reproduces observed trends in the compositional 
variation of o'(x), including the observation of a step- 
like jump as glasses enter the Intermediate phase near 
Xc(l)=9%, and an exponential increase when glasses be- 
come flexible at x > Xc(2)=38%. Since is found to 
be small compared to network strain energy (E^), we 
conclude that free carrier concentrations increase linearly 
with nominal Agl concentrations, while the super-linear 
variation of fast-ion conduction can be traced to changes 
in carrier mobility induced by an elastic softening of net- 
work structure. By better understanding the role played 
by the " alloyed network" on fast-ion conduction in a spe- 
ciflc Ag based solid electrolyte glass, we have also taken a 
step towards better understanding ion-transport in sev- 
eral other types of solid electrolyte systems, including 
modified oxides and modified chalcogenides, systems that 
are also expected to display a parallel conductivity re- 
sponse as a function of network connectivity. Variations 
in /?(x) reveal a square- well like variation with walls lo- 
calized near Xc(l) and Xc(2), and suggest filamentary ( 
quasi ID) conduction in the Intermediate Phase ( Xc(l) 
< X < Xc(2)) but a dimensionality larger than 1 outside 
the IP. The present findings, characteristic of dry sam- 
ples and showing two elastic thresholds, may possibly go 
over to only one elastic threshold in wet samples in which 
the IP completely collapses as suggested by earlier work 
in the field. 
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